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ABSTRACT 



This report investigates subsonic diffuser performance, with 
emphasis on conical and annular geometries. A correlation is 
presented which aids in the prediction of performance. Two an- 
nular diffusers were designed and tested to help substantiate the 
correlation . 
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I. INTRODUCTION 



The diffuser is a device used in fluid mechanical systems to 
convert kinetic energy into static pressure by decelerating the 
flow. There are several basic straight-walled diffuser shapes 
such as conical, two dimensional and annular as shown in Figure 1. 
Much attention has been given to two dimensional diffusers and 
several successful correlations have been developed, but annular 
and conical diffusers have not been as exhaustively studied and 
good new correlations would be very useful to the designer. 

A. A GENERAL DESCRIPTION OF THE FLOW 

A diffuser should convert kinetic energy of the flow into 
pressure with a minimum of losses due to viscous effects. Also it 
is desirable to make optimum use of the geometric area ratio, 
because the area ratio prescribes the amount of diffusion that is 
possible under ideal conditions. Bernoulli’s equation relates 
the total fluid energy to the static pressure and kinetic energy; 
and if the flow experiences no losses due to friction or other 
non isentropic processes the total energy of the fluid remains 
constant. Any decrease in velocity will cause the static pressure 
to rise, so that in an ideal diffuser all kinetic energy losses 
show up as pressure gains. Unfortunately, most fluids experience 
losses, and the goal of any flow analysis is to keep track of the 
fluid energy whether it is dissipated or converted from flow 
energy into static pressure. Dissipation of kinetic energy takes 
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place inside the boundary layers where large velocity gradients 
exist so that the viscous effects become important. 



The boundary layer is a thin region adjacent to the wall in 
which the velocity increases from zero to the free stream value. 
Since the fluid in the boundary layer is moving slower than that 
in the free stream, it does not have as much kinetic energy as the 
bulk of the flow. The transverse velocity gradients that are 
present give rise to shear stresses which are related to the veloc- 
ity gradients through the coefficient of viscosity for laminar 
flow. In a turbulent boundary layer the shear stresses are not 
related to the gradients in such a simple fashion. However, it is 
possible to model turbulent flow as an average motion plus time 
dependent fluctuations whose average value is zero. When the 
boundary layer equations arc modified for turbulent flows , addi- 



tional shear stress terms are present known as Reynolds stresses. 
These Reynolds stresses are in fact the predominant stresses for a 
turbulent boundary layer. 

Since the flow in diffusers is constantly decelerating, due to 
increasing area, the flow faces an increasing pressure as it 
traverses the channel. To overcome the increasing pressure the 
flow must transform kinetic energy into pressure. However, as 
mentioned the flow in the boundary layer is not as energetic as 
the free stream, but it feels the same longitudinal pressure gradi- 
ent. When the fluid near the wall does not have the energy to 
overcome the pressure rise, the flow higher up in the boundary 
layer must pull it along by shearing action. In this sense the 
turbulent boundary layer is more capable of meeting an adverse 
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pressure gradient than a laminar one because of the additional 
Reynolds stresses. For this reason it is advantageous for dif- 
fusing flows to be turbulent in order to withstand the decelera- 
tions imposed in a diffuser. A point is reached, however; where 
the outer layers can no longer pull the inner layers, and the flow 
is forced away from the wall and separation results. Large scale 
separation with recirculating flow causes intolerable losses in a 
diffuser and is to be avoided. Any time the boundary layer does 
not have sufficient energy to overcome the pressure gradient it 
will separate, so it is necessary to find out something about the 
way the pressure gradient behaves in a diverging channel. 
Bernoulli's equation is written as follows 

P t = P + i V a) 

Since the boundary layer feels the same pressure distribution as 
the free stream, it is only necessary to calculate the free stream 
pressure distribution which shall be calculated here assuming 
inviscid one dimensional flow in the core. Taking the derivitive 
of (1) gives 

dp/c/\ r ~ ydv/d\ (*) 

with continuity and constant density 

m = -fAV = -f, u) 

or 

d A/a = v ’/= (- v7a 'W A co 



11 



replacing (4) in (2) 



dp/Jx = (-fVVA) 



( 5 ) 



with continuity 



Jp/Jx = (fd* V, 1 / A 3 ) JA/Jx 



( 6 ) 



or 



dWifrt/dx = dO/$)/d* * (z4,Va 3 ) c/A/dx (?) 



(7) shows that the pressure gradient is the strongest at the inlet 



constant for a given straight walled diffuser. This is compatible 
with the state of the boundary layer which is "young” at the en- 
trance of the diffuser and is able to overcome large pressure gra- 
dients. The large initial gradient suggests that a boundary layer 
should be as thin as possible at the entrance to a diffuser for 
best performance; i.e., avoidance of early separation. This fact 
has been substantiated by many researchers by increasing the inlet 
boundary layer thickness and noting a decrease in performance, 

Kline ^Ref. J- Returning to the expression for the pressure gra- 
dient, (7), it can be seen that the gradient may be high initially. 
But a large area gradient also insures that the pressure gradient 
term will decrease rapidly, because the denominator of (7) contains 
the local area raised to the third power. Although a general 
description of the flow conditions in a diffuser has been given, 
the task remains to determine losses due to the factors considered 
above . 



of a diffuser since the area 



gradient, , is very nearly 
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B. REVIEW OF PREVIOUS WORK 



Many attempts over the years have been made to explain and 
predict the actions of the flow in the diffuser. As the previous 
section suggests, it is not a simple process. The major stumbling 
block is that the phenomenon of turbulence is not well understood, 
and to predict separation the growth of the boundary layer must be 
calculated. In order to bypass the turbulent flow calculations, 
more empirical methods and simplifications have been used to pre- 
dict performance. One of the simpler techniques is that of the 
equivalent cone angle method as described in Gleason [Ref. {]. 
Basically, all flow cross sections are related to an equivalent 
conical flow with a corresponding wall divergence angle. If the 
wall divergence angle, which is a direct function of pressure gra- 
dient, is below a certain value the diffuser is judged sound. 

This method is widely used for design but has not always proved 
adequate. The next level of effort has been directed towards gen- 
erating performance plots. 

After testing a sufficient number of diffusers of a given type, 
constant pressure recovery lines were plotted on a graph of perti- 
nent geometric parameters such as wall angle versus length. 

Kline ftef. Q has done extensive work in this area for two dimen- 
sional diffusers and Sovran and Klomp [Ref. have done the same 
thing for annular diffusers. Through flow visualization techniques 



flow regimes for different two dimensional geometries. Figure 6. 
The line of first stall, a-a, is presented as the optimum perform- 
ance line, because at or just after this point the boundary layer 



Kline 




showing the expected 
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has overcome as much of the pressure gradient as possible and is 
nearing large scale separation or transitory stall. No similar 
flow regime graphs have been developed for the annular or conical 
geometries 

Several years ago Kline ^Ref. l] developed a more analytic 
approach to the problem. He has attempted to solve the turbulent 
equations in predicting performance. The boundary layer and con- 
tinuity equations have been combined into a mementum integral 
equation and used with several shear correlations to work up a set 
of five linear first order differential equations with non linear 
coefficients. However, the solutions of these equations alone do 
not predict the performance, a stall criterion is needed. Kline 
has combined several parameters including boundary layer thickness 
and rate of area growth into a parameter. When this parameter 
decreases to a certain value stall is said to occur. The param- 
eter has no physical basis, but it does correlate line a-a of the 
flow regime chart. In Reference 1 many reports have been analyzed 
and the criterion in fairly successful for 2D-diffusers . The 
variation in performance with inlet boundary layer thickness -is 
also demonstrated. Though Kline appears to be fairly successful, 
the method is rather approximative and requires the solution of a 
large set of equations. However, Vavra Qlef. 8j[ has developed a 
performance parameter which is easy to calculate from geometry and 
can be derived from basic principles. 
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CONICAL 



DIFFUSER GEOMETRIES 




FIGURE I 




FIGURE I (CONCLUDED) 






II. DERIVATION OF THE SHAPE PARAMETER OMEGA 



Vavra [Ref. Sj presented the basic derivation, and Gapp ^Ref . 3^ 
did a thorough development which is presented below. 



The momentum theorum 



f div\ s V 1 - f (Iwjr V, - - f * f dj- { it ds 

d) * - (b • s s 



( 8 ) 



was applied to the differential flow element of Figure 2 to yield 






p/I - fp+dp) (A4dA)~H 3 p<IA 

» ^ ^ ^ 

After expanding and rearranging terms where «r» ( and “ ^ are the 

outward pointing normals of each respective side and •L is the 
tangent vector, the theorum gives 



^ AV dtf « - A^f~ 2*C dl 



rearranging 



dp = - -f '/ dY' - 2 * C d L/A 



( 10 ) 



(ii) 



or 



d ( V 1 / + dp/-/ - C - 2'/f)( c /A) dU 



( 12 ) 



Equation (12) may be derived for an annular diffuser as well. 
From Figure 2 (a) for an annulus 

cf L, = d L c. oz xs ( l3 ) 

A, = A/Cos. V (>’> 
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( 15 ) 



From Figure 2 (b) 

c = nr ( 5+ rj = H-rr r„ 

Also c'=<M- r„ for small divergences. Replacing (13), (14), and 
(IS) in (12) 

d ( 1 / 72 .) +■ d - ( - Zr/-f) ( c'/A , cos is)<SL l CM>f 



T C dL, /-f A , 



So finally considering only sections normal to the axis and 
omitting the prime 



d ( V V 2 -) + <*?>// 3 (-?/■/) <C/A) du 



as before. From the first law of thermodynamics 



J r\ 
U 



J 

u 






r 



u 



xr rr f. 



1 / 



. ^ t/ ~ I . _ 

»/ a p - ok- y< a p - 

f / 1 



nr - 

* ^ 



U-Oj 



or 



dp/j = cIL-Td* 



Replacing (17) in (12) gives 

dL+ J ( v 1 /?.) - T ds » (- 2-/f ) C c/*) <«- 



( 17 ) 



(18) 



combining the enthalpy and velocity derivatives into the total 
enthalpy form 

d+l = T<f.s - C ?/f) (c/A) dL (19) 

Since the idealized one dimensional process is adiabatic, with no 
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energy input through the walls, the total enthalpy is constant and 
its derivative is zero. Therefore (19) reduces to 



Tds = (t/f) (c/a) AL 



( 20 ) 



The shear losses in a turbulent boundary layer are not easily 
determined, but they do give an indication of overall losses. An 
average shear stress coefficient was defined for the entire chan- 
nel to reflect these losses. 



C jf. - Z'M/l-fV- 

or using continuity 

” (Qf Z 2 *) **'» "/-/ A*“ 

Replacing (22) in (21) yields 

Tds ~ (C^/ 2 .) ) C6L 



( 21 ) 



( 22 ) 



(23) 



Now a similar differential entropy change will be developed from 
purely thermodynamic considerations for an adiabatic poly tropic 
compression. Referring to Figure 3, the polytropic efficiency is 
defined to be 7^= 41: j;or using the isentropic relations for a 



dT 

differential compression 



(t •+ dT u )/T = ( P+ dp/p) 






(24) 



Expanding in a binomial series and neglecting higher order terms 

V/£-f 

(l *V - I + .Jl 4 TZr-t«s ( 25 ) 



Replacing (25) in (24) and rearranging 






dT,s / dr = (r/dr)(d ? /p)if 



(26) 



19 



using 7 p - dr,,/c/T 

iT/T - (k-l) /&/> X) d?/ P ( 27 J 

Introducing the perfect gas law and Cp - ^ nto t * 1G first law 

of thermodynamics yields - 

Tds = R ( V/;r-0 JT - (n.T/p) dp 

or 

d.% - R Q ( X/X-i) ciT/r ~ ( 28 ) 

The final differential entropy form results by using ( 27 ) 
in ( 28 ) 

ds/R ~ ( V'fp -l) dp /p ( 29 > 

In order to combine the geometric and thermodynamic developments 
the differential entropy changes. ( 29 ) and ( 23 ). can be equated 

(Cr/zr) .tid- cdL ~ RC'yp-t) dp/p W 

Introducing the dimensionless referred mass flow rate 

into ( 3 CT) yields ' 

Cc^/zr) ( M,.'“ A, 5 - p, 1 / £T , I s - a- ) CdL " 



(T / / ) ( V^tp - 1 ) dp / p 

again using the perfect gas law -f = FY lt 

(?,/j)(Q /Z.T) A }/A l \ cdL = 

Cr/f) ( '.'-np -i) dp, 



( 31 ) 



( 32 ) 
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rearranging 



C^/2.) ( h * r~ CdL = (f/f,)(fy-M(*/P,) 



(33) 



For a polytropic process fyj* 1 = constant: so (33) becomes 

(<#/*.) * r >- ( A.VA*) c d L - ( l /"lf -i)frp, )' l 'd P'n ) (34) 

The equation is now in a suitable form to be integrated. Vavra 
has called the left side ax, and integrated it from the channel 
entrance to exit along s streamline so that the constant total 



enthalpy assumption is valid. 

L 



J 2x, = J ( c - ( a* /A 1 ) CJL 



(35) 



As previously . discussed is a constant representing, shear losses 
and may be brought outside the integral. Also referred flow is a 
constant for a given flow rate . 

r ju « / r , \ . * 1- r 'A ~ * » 3 » j % 

J aX, - ( C -j / 2L. ) i K * fn-yy-ou 



ww 



The integral of (36) contains geometric characteristics of the 
diffuser and has been defined as the shape parameter omega (*Xw) . 

SI - f (At/AY dL ( 37 ) 



The product G, must be related to diffuser performance for 

T 

SL to be useful. In order to do this the right side of (34) 
must also be integrated. 

dK - ( '/>-/) C */?,) Vk dC p/p,) 



Integrating as before from entrance to exit. 

X. = f ^ ^ - 1 ) ( P/P >,)* d(tyr.) 



(38) 
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( 39 ) 



By assuming constant polytropic efficiency and n constant 

= (fyp -l) (h/ n *i) ( ( fyp,) V ~0 

For the differential change |> r p ( 4 . p with &p 4 <£p f and 

again using a binomial expansion 

n+ i 

[ (?,+ 4p)/p,] n - I + 0 +W + higher order terms 

replacing in (39) 

K i. - ( V'n? “0 ^p/p. tw) 

equating (36) and (40) 

£(a-/a) s ( c m,)cIl ~ 

■( c */i)v" r ' k JL = <f '/^p "0 6P/P, w 

substituting for according Lu Lhe definition 

d, V,.{£t?/ P, a, 

gives 

(C + /2.) fp I/, 1 P-T t ) -ft /<V A,'1 =(yj/>~l) 4P/P, 

or 

iSbC 4 - C IA\>/ i, V', 7 -) c '/-if -1) f42) 

Since constant total enthalpy was assumed there is no change in 
total temperature so rewriting (28) 

is 3 - R. ( of Pt /Pt) 

equating (29) and (28a) 

- d P-t / Bt - ( '/-v 



(28a) 

(29a) 
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(29a) is in a suitable form to be integrated and is a straight- 



forward log. 



(-U ^ -0 l I - *1. 






or 



In ( P-fc, /h ) - ( Vyp -t ) 



(43) 



again using 



h, /P ±L = and r ,+ 4 

with the usual expansion for logs; i.e. ln(l+x) = x, (43) becomes 



4 ?/P, ( /P-ty 



(44) 



or 



(■ /2p = - ^p t ( ~ & §./0+ $ M I 1 ) t" 5 ) 

All flow considered was incompressible with so 

C '/'np -i) df> - -a 

replacing (46) in (42) 

-n. Cf » - i a p + r p* - POAt 7; kV 



(46) 



(47) 



or 






A/C. 



(48) 
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( 49 ) 



However Cp = coefficient of pressure 

c p * C 

c n = pressure recovery if the process followed the P, to 

_ line of Figure 3', i.e. ideal. 

* isen 

c ?i = /(tfvr) (so) 

For isentropic flow = O so (50) becomes 

Cp. = ( -//i'/, *• •• /A = I - ( A '/I\S (SOa) 

Replacing (50a) in (37) gives 

(P* = 

[ P, + i / \Z, ~ i -f K ^ ( i - Cp-)3 /(£ / V, (si) 

using (49) and (48) in (51) yields 

SiCr. ~ C ?i.“Cp (52) 

(52) is the basic result of the theory. The amount a flow 
departs from isentropic conditions has been shown to be a product 
of the geometric parameter omega and an average shear stress 
coefficient . 
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FIGURE 2 (c ) 

ANNULAR CHANNEL SECTIONS 




FIGURE 2 (b) 



26 




<^aa 3 

FIGURE 3 



f-s diagram foran adiabatic comp re ss ioh 



III. ORIGINAL SUGGESTION FOR USE OF OMEGA 



Once equation (52) was derived, Vavra {^Ref. 8^J investigated 
some two dimensional data by plotting the losses as a function of 
the shape parameter. His plot is reproduced in Figure 4. It can 
be seen that the losses tend towards a minimum value in the neigh- 
borhood of an omega of ten. For this reason Vavra suggested that 
an omega of ten might be a good diffuser design value. 

In an earlier investigation of the shape parameter, 

Gapp £"kef. 3j[ reduced the* annular data of Sovran and Klomp f Ref . bj 
so that a plot similar to Figure 4 could be generated as shown in 
Figure 5. Although the plot shows the same trends as Vavra 's 
graph, there does not seem to be a value of omega to insure mini- 
mum losses. Additional two dimensional data of Kline Q\ef. s7 was 
also plotted and has very similar characteristics, but does not 
offer any firm evidence that one value of omega is to be strongly 
preferred over any other. A somewhat different approach was to 
plot lines of constant omega on the two dimensional flow regime 
chart Fox and Kline ^ Ref . 5*], Figure 6. It was seen, however, that 
a given value of omega did not necessarily indicate the state of 
the flow. Although the above evidence indicates that Vavra’ s 
original suggestion of how best to represent losses as a function 
of omega was not the best implementation of the shape parameter, it 
does not adversely reflect on the validity of omega as a good 
parameter. Omega has a definite physical significance which is 
fundamentally different from previously suggested parameters 
embodying all pertinent geometric characteristics. A der. .'utim 
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of the physical meaning of omega has already been given and an 
additional aspect of omega will be given in a later section. 

The effort to achieve small losses as reflected in a minimum 
value of JIC+ may not always be the chief goal of a particular 
diffuser. The ideal pressure recovery coefficient is a function of 
area ratio and as such is a measure of the maximum diffusion ob- 
tainable for a particulr geometry. If the designer wishes a dif- 
fuser to convert a large portion of the flow evergy into pressure, 
a larger area ratio is called for. However if only a modest rise 
is called for with a fairly uniform exit profile the reverse is 
true . The application calling for a large pressure rise may be 
able to tolerate higher losses when maximizing recovery while the 
need for a uniform velocity profile will tend to require minimum 
losses In proposing a theory or correlation for performance the 
two above conditions must be considered so that an optimum design 
criterion may be formulated. Before investigating diffuser 
characteristics further a somewhat simpler derivation of omega will 
be given to relate the shape parameter to diffuser performance in 
a different way than previously derived. 
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F I G UR E 



IV. AN ALTERNATE DERIVATION OF OMEGA 



The non dimensional group used to define omega was a logical 
consequence of the combination of the momentum theorem for an 
arbitrary channel and the expression for a differential entropy 
rise in an adiabatic poly tropic compression. The derivation was 
necessarily lengthy and the expression for omega may be obtained 
from a more basic approach. 

For the same one dimensional channel flow of Figure 2 



Taking the derivative of (53) with respect to length and assuming 
constant total pressure 




(53) 




(54) 




o 



(55) 



and from continuity for constant density 



vda = - AiV 



or 



V cl \J - C “ V*~/A ) d A 



(56) 



replacing (56) in (55) gives 




( 57 ) 



and using the definition of mass flow from (54) 



rearranging 

< ip /({ fV, x ) = 2- C a,/a)* fM,) cm/jl.) jl ( 59 ) 

and integrating both sides 

C]>^ ~ f C^i/AY 0/A,) ( d A/dl) Jl m ( 60 ) 

The integral in (60) looks very much like the definition of omega, 
but the area rate of change must be related to the wetted perimeter 
for each diffuser type. 

A. CONICAL 



A Oh = 


2 rr 


^ R-, * 


L TA N 0) 2, 


(61) 


dA/dL — 


2 


rr tA-V 


0 ( pL LTAa/O) 


(62) 


C(l.) = 


2. fj~ ( fiL t + 


L T A A/ &) 


(63) 


combining (63) 


and 


(62) yields 






= 


(d 


A / d L) 


/ T A sj O 


(64) 


and replacing 


(64) 


in (60) 







Cj>: s 2. f C(l) t*«/ 0 <iL -XTti'VQiJlAte) 

1 a 



Therefore for the conical case the ideal pressure rise is a 
function of the wall angle and omega. The factor relating 
Cp ^ to omega has been defined as PAR4C . 
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B. TWO DIMENSIONAL 



M/dL « 2&,TA A/© ( 66 > 

C^U) a 2-6, -t- 2W, + 4 LT Af/o - 2.1, + 2 WCL) 

or 

2B, s CO.) -2 IVfL) (67) 

combining (66) and (67) 

d M/c/ L - C( l) V# ( I ~ Z. W( L ) / Cc'l) ) (68) 



^Pi ~ / 2 - r ^ A/ 0 ^ / ~ 2~*V(l\/C{l))(~±) (^~)u'L(69) 

Equation (69) is also similar to the expression for omega, but an 
investigation has shown that Cp^ is not related to omega by a 
simple constant analogous to that of the conical diffuser. 

C . ANNULAR 

An annulus is a complicated shape, but it shares certain 
fundamental geometric characteristics of the cone such as axial 
symmetry. In fact, an annular diffuser is nothing more than a 
conical diffuser with an axisymmetric center body. For an annular 
diffuser: 

cM/o'L - 2tfL (TA d/- & 0 -TAa/*” &i) 

zrr ( an rrfv&o - fz tf-l TAdOi) (70) 

CL Cl) - ITT C (71) 
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dt\/d L "" f( T^a/ B 6 Rt,(L) -TAa/q* fSjc ^£-))l C(l) (72) 
L* * + r-o O.) ) 



By examining the factor relating area gradient to the local cir- 
cumference in (72) , it can be seen that an expression similar to 
(65) can be derived for annular diffusers. The factor is defined 
as PAR4A . 



PAfc-HA^ 12oM -TA*/&i P-rCOl / Qi'.rO.) 



By referring to Figure 7 and making use of the theorum relating 
similar triangles from geometry it can be seen that PAR4A may also 
be expressed as a function of the inlet geometry. 



RRAT = inlet radius ratio = l ■ I 

It can be seen that PAR4A is constant for a given annular diffuser 
family. A family is defined as a set of diffusers with constant inner 
and outer wall angles and the same radius ratio. For example all 
diffusers with wall angles of 15 and 20 degrees with a radius 
ratio of one half would comprise one family while all diffusers 
with wall angles of 10 and 25 degrees with a radius ratio of .7 
would comprise another family. The previous derivation may be 
used in (52) 




/ £>H > / L,' \ L 



~ (TAaJ& 0 ~£.fb/4r T&'V&i) /( /•#“ 



-fi- = ^ pC ~ ^ P 



(52) 
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or 



UhCj. ~ ( i - Cp/Cpi) 

P/tpc) ~ 




(52a) 



Although the previous derivation of omega does not allow direct 
calculation of performance, it does give a better understanding of 
the shape parameter. Once the above derivations were formulated 
a parameter study was undertaken to observe the response of omega 
to certain geometric factors for each diffuser shape. 
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FIGURE 7 
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V. PARAMETER STUDY 



The parameter study was done with the aid of several computer 
programs which are listed .and explained in Appendix A. To do a 
complete study would have been a major undertaking with all the 
associated variables of conical, two dimensional, and annular dif- 
fusers. However, a representative study was done by varying each 
parameter a reasonable number of times to establish trends. The 
conical diffuser was the simplest shape and was studied first. 

The conical variables were taken as inlet area, length, and 
wall angle; but to have a manageable amount of data the inlet area 
was set equal to unity and the other parameters were varied. 

Figure 8 shows that omega is a strong function of wall angle and 
past a certain non dimensional length, the wall angle dominates. 
For a high divergence angle the characteristic length is short be- 
fore a constant omega is achieved while the opposite is true for 
the lower divergence angles. The wall angle is inversely propor- 
tional to the final value of omega as shown in Figure 8. An 
interesting fact about Figure 8 is that an omega of ten results 
for a wall angle of about three degrees. Conical diffusers with 
small divergence angles have low losses, but they do not neces- 
sarily have the highest pressure recovery for lengths common in 
practice. A better look at omega in relation to flow conditions 
was found by plotting constant omega lines on the two dimensional 



Figure 6 was briefly mentioned earlier, but there arc several 
aspects of the plot that were not discussed. For the two 



flow regime chart of Fox and Kline 




39 



dimensional geometry, omega was found to be dependent on angle, 
length, and inlet geometry. As Figure 6 demonstrates, a particu- 
lar value of omega in the two dimensional regime can indicate a 
full range of flow conditions from unstalled to large transitory 
stall. In addition to not indicating conditions in the two 
dimensional diffuser, the constant omega line takes on different 
values depending upon the inlet geometry. 

Finally the characteristics of annular diffusers were studied. 
The annular diffuser is the shape of most interest for turbo- 
machinery applications and unfortunately is the least understood. 
There are numerous parameters to vary: wall length, radius dif- 

ference, mean radius, inner wall angle, and outer wall angle. 

Due to the large number of variables no plots analogous to 
Figure 6 and Figure 7 were apparent to show omega trends. However, 
one graph which combined inlet geometry and wall angles was plot- 
ted versus non dimensional length. Figure 9 has the advantage of 
incorporating as many geometric parameters as possible into a 
meaningful group by plotting area rate of change versus length at 
a constant value of omega. The area rate of change is about con- 
stant for a given set of wall angles and inlet geometry as can be 
seen from (70) when ( T A*/ 1 & 6 -T /^ /V^^bjean be neglected, 
and it is also related to the pressure gradient. The plot shows 
the same general trends as the other two geometries by indicating 
the omega dependence on area rate of change and length. The 
figures all show that high omegas can only be achieved with low 
area rate of change, but several geometric conditions may result 
in low values. A diffuser may be very short or have a hi h area 
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gradient in the case of a low omega value. Although the parameter 
study indicated geometric trends, it did not offer any insights 
into performance. To gain an understanding of diffuser perform- 
ance the published data was analyzed more closely. 
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VI . DATA CORRELATIONS 

Previous attempts to predict or explain diffuser performance 
relied on a number of parameters and graphs that were rather 
unweildy to use, or an oversimplified approach such as the equiv- 
alent cone angle method was employed. Gapp ("kef . 3^ looked at the 
data of Sovran and Klomp [Rof. 6j| in the manner suggested by 
Vavra (kef. and tested several diffusers, but he was unable to 
develop a method of prediction based on his work. It has been 
shown in previous sections that the parameter omega has physical 
significance in relation to diffusing flows, and it was reasonable 
to assume that omega is in some way related to performance. If 
some means of determining the shear term, C^, could be found; 
equation (52) would give performance directly. Several attempts 
were made to arrive at a technique for determining C^, but some 
other methods of predicting performance from curves bear mention- 
ing first. 

If ~T\. is small the recovery is very close to ideal, but 

must be minimized for a given area ratio to insure maximum 
pressure recovery. This fact suggested a possible utilization of 
the original plot of Vavra . By drawing lines of constant area 
ratio on an versus Cf graph, an optimum omega range can 

be determined at the given area ratio. For the given area ratio, 
the omega range would guarantee minimum losses and maximum pres- 
sure recovery close to the desired one. A plot such as Figure 10, 
which is a line of area I'atio two from the data of Sovran and 
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Klomp [Ref. bj, would result. It shows an omega range of five to 
eight for minimum losses. Unfortunately Sovran and Klomp did not 
have a wide range of area ratios, and to do the required testing 
to generate other data would be very costly and time consuming. 

Another plot was discovered, but it did not have a wide 
range of applicability. However, it is worth mentioning since it 
correlated a large amount of the annular data of Sovran and Klomp. 
The non dimensional area gradient, ( c/$/cf L ) /A P* , is related to 
the pressure gradient, as shown in (S7) . After several attempts 
to relate the area gradient to losses a graph such as Figure 11 
was plotted. It is a plot of non dimensional area gradient versus 
losses for a given non dimensional length, and a straight line 
resulted for a range of lengths . The resultant straight lines 
suggested plots such us Figures 12 unci 13 in un ^ittsrnpt to 
relate all the annular data. Figure 12(a) is a list of symbols 
used in Figures 12 - 15. However, considerable scatter is evident 
in addition to some definite trends'. Two distinct lines are vis- 
ible and some data considerably off both curves. The two distinct 
lines correspond to families with inlet radius ratios of .55 and 
.70 respectively. The data below the lines are from families with 
inner and outer wall angles that' differ by no more than two degrees 
These plots are of limited value at best and are rendered even 
more questionable because the diffusers with the best recovery, 
the families with similar inner and outer wall angles, do not even 
fall on the curves . 
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The independent derivation of omega introduced a geometric 
parameter that related the ideal recovery to omega (see 
equations (52) through (72) . The parameter, PAR4, had a different 
form for each class of diffuser, but in all cases it corresponded 
to the quotient of the area gradient and circumference for a 
given axial station. If the shear stress coefficient could be 
determined, the theory of Vavra would give performance directly 
from equation (52) . To see if the parameter, which has been desig- 
nated PAR4A and PAR l JC, had any relation to the shear stress coef- 
ficient, Figure 14 was plotted with the conical data of Cockrell 
and Markland C Ref - 2 J and annualar data from Sovran and Klomp. 

The conical data consisted of twenty-four diffusers which were 
tested at two boundary layer thicknesses and approximately eighty 
annular diffusers. No information was given concerning the exact 
dimensions of the inlet boundary layer for the annular tests, but 
all the annular diffusers were tested on the same apparatus with a 
thin inlet boundary layer. The set of conical diffusers which is 
displaced upwards on the curve had a thicker inlet boundary layer 
and exhibited the expected higher losses. Some scatter is evident 
in the figure for the longer annular diffusers, but a definite 
curve is discernable . 

Once Figure 14 was established, the next logical step was 
Figure 15. To enable the direct determination of recovery from 
a known omega, PAR4 was plotted versus the quotient of pressure 
recovery and omega. Again a good curve with small scatter resulted. 
The conical data with the thick inlet boundary layers were dis- 
placed downwards as expected. Figure 15 differs from the previous 
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annular correlation in that the data fall close to the line, and 
the wide scatter attributed to families with similar inner and 
outer wall angles is not present. Also the curve does not seem to 
depend on inlet radius ratio. With Figure 15 and a knowledge of 
geometry for a conical or annular diffuser with a thin inlet 
boundary layer thickness, the pressure recovery can be determined. 

Besides predicting the recovery, the correlation presented in 
Figure 14 can be utilized in the development of optimization 
criteria. The designer is often faced with the problem of maxi- 
mizing the recovery for a given .length. The annular case will be 
discussed, but the development for a conical diffuser is similar. 
It was shown that 

SLC ^ ~ Cpj “ C p and C = P tl. 4 SL> 



Combining the equations gives 

Cp = Cj>- ( I - Cx. /Mp-H) ( 73 ) 

However, Figure 14 shows that is a function of PAR4 . The curve 

is of the form C r ** , where ci. and b are 

V 

constants that can be determined from polynomial curve fitting 
techniques used in numerical analysis. Also the ideal recovery 
coefficient can be shown to be a function of inlet conditions and 
non dimensional length. 

Cp. * I - fa/Aj*- 




C Hr I -h P-H l) / C & 77 * flN/ + * £ ( L /& K 

( 7V?// z 0 o “ Tsirf’Q i ) f ' :C,r ' ( " ' * 

fc HI TA<\t 0 ' u 1 



(74) 
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By replacing (74) in the expression for the ideal recovery coef- 
ficient and combining this result with an expression for Cf as 
suggested, the pressure coefficient will be expressed as a function 
of inlet geometry and non dimensional length. When an expression 
has been developed, the theory of maxima and minima of a function 
may be used. By taking the derivative of Cp, with respect to non 
dimensional length and setting the resultant expression equal to 
zero,- the recovery could be maximized. To obtain a given recovery 
in a minimum length, the correlation could be employed in a similar 
fashion . 
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VII. EXPERIMENTAL WORK 



A. EXPERIMENTAL EQUIPMENT 

Two annular diffusers- were designed and tested in order to 
provide further data in the evaluation of any trends or correla- 
tions discovered from analysis of the literature. The diffusers 
were designed with the aid of program DESIGN which is listed in 
the appendix. With a given inner wall angle, fifteen degrees, the 
program varied the outer wall angle and length to maintain a given 
area ratio and calculated omega for each case. Both diffusers had 
an area ratio of three and an inner wall angle of fifteen degrees. 

Model one had an axial length of 7.73 inches with an outer wall 

angle of nineteen and one half degrees, and model two had an outer 
wall angle of eighteen degrees wiLh an axial length of 0.36 inches 
Both diffusers had an inlet hub radius, Rhl, of 1.875 inches and 
an outer inlet radius, Rtl, of 3.123 inches. Model one had an 
omega of four and model two had an omega of five. Both were con- 
structed from phenolic resin. The above information is summarized 

in Figure 16. Figure 17 is a photograph of the inner body and 
model two. In order to supply a uniform flow to the diffusers, a 
contraction cone was designed by the method described in Ref. 7 
and was also fabricated from phenolic resin. Figure 18 is a photo 
graph of the contraction cone and model one. The inner body was 
supported by three struts mounted in the cone as can be seen from 
the overall drawing of the assembly in Figure. 19. 
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Static pressure taps were placed at intervals of two inches 
with each station having three taps one hundred and twenty degrees 
apart on the inner and outer bodies. Two stations one hundred and 
twenty degrees apart were also provided one half inch from the in- 
let plane to made radial surveys. A survey station was also lo- 
cated 8.25 inches upstream from the diffuser inlet station to make 
cross sectional velocity provile measurements and detailed measure 
ments of the boundary layer. Hot wire anemometer techniques were 
used for these measurements and the static pressures were read 
from a manometer board. A United Sensor cobra probe was used to 
measure the inlet dynamic head and it is shown with the hot wire 
and model in an overall view of the assembly in Figure 20. 

The flow delivery system consisted of an axial compressor, two 
settling chambers, flow straighteners , and an eight inch pipe. 

The compressor was an Allis-Chalmers twelve stage axial compres- 
sor which operated at a pressure ratio of three to one. After 
leaving the compressor, the air was' cooled and sent through flow 
straighteners into a large plenum. After passing through a flow 
straightener the air passed through a sharp edgcd-orifice for flow 
measurements. Before being supplied to the eight inch pipe, the 
air passed through the second plenum. After leaving the eight 
inch pipe the flow entered the contraction cone and then the model 
Meriam micromanometers indicated the pressure readings for flow 
rate calculations. 
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B. EXPERIMENTAL RESULTS 

Any irregularities in the inlet velocity progile to a diffuser 
can cause separation and poor performance, therefore it was neces- 
sary to insure that the inlet profile was uniform. The hot wire 
was used to make a survey at the measuring station on the contrac- 
tion cone to measure the profile and boundary layer thickness. 

The contraction device provided a uniform inlet flow, as can be 
seen by the measurements of the velocity distributions at the sta- 
tion 8.25 inches upstream of the diffuser inlet presented in 
Figure 21. The boundary layer thickness was about .2 inches and 
the displacement thickness was calculated to be .051 inches and 
the momentum thickness .014 inches. The flow at the inlet plane 
of the diffuser is also important. It was checked with a hot wire 
and the velocity profiles for both models we r> e found hr> he rather 
flat as shov/n in Figure 22 . The assurance of uniform inlet flow 
was important for the possibility of comparing this data with the 
other existing data. Uncontrolled non-uniformities resulted in 
flow separation at the diffuser inlet in previous tests with a 
badly designed contraction device. 

The most important result of a diffuser test is the pressure 
recovery of the models. The static taps allowed a measurement of 
pressures around the inner and outer bodies and the cobra probe 
was used to measure the inlet dynamic head. The pressure coef- 
ficient is defined as 

c p --- O /%, 

where K is taken as atmospheric pressure. Both diffusers had an 

& 

ideal recovery doefficient of .989. The non-dimensional sccitic 
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pressure distribution for model one is shown in Figure 23 . The 
recovery and losses are also shown on the graph. Model two also 
had good recovery as shown in Figure 24. The tests are presented 
in Figures .14 and 15. It can be seen that they agree with the 
published data presented in these figures. 
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Figure 17 Inner body ond model two 




p^gure 18 Contraction cone and model one 
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Photograph of equipment 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



Diffuser performance was investigated in order to obtain a 
simple correlation to use for design purposes which related dif- 
fuser geometry to performance. From analysis of the available 
experimental data a geometric parameter suggested by Vavra, 

-A- (Section II, p. 17), was used to relate the important per- 
formance quantities to geometry in a simple manner. A correlation 
relating the available experimental conical and annular diffuser 
data was finally established between diffuser geometry and per- 
formance without, considering the status of the inlet boundary 
layer. The method of utilization of the correlation in minimizing 
the diffuser length for a specified pressure recovery, and maxi- 
mizing recovery for an available length, was also indicated. 

Two diffusers were designed and tested to provide additional 
experimental data and results were found to agree with the 
developed correlation. More diffusers should be designed and 
tested in the region of the correlation where the data are not 
plentiful. The influence of the inlet boundary layer and Reynolds 
number must be investigated, since it is felt that these two fac- 
tors are among the causes of the scatter in the derived correlation. 
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APPENDIX A 



Program Design was formulated to aid in the design of the 
diffusers mentioned in the experimental section. The purpose of 
the program was to vary the outer wall angle from sixteen to 
twenty-two degrees in increments of one quarter while always main- 
taining a specified area ratio. This procedure was done for 
three different area ratios. The program calculated omega for 
each wall angle combination by the use of a numerical integration 
scheme which employed an external function. 

List of Symbols 

A1 - inlet area 

A2 - outlet area 

A - area 

C - circumference 

I - logical control variable 

RH1 - inlet inner radius 

RT1 - inlet outer radius 

AHR - inner wall angle in radians 

ATR - outer wall angle in radians 

RT2 - outlet outer radius 

RH2 - outlet inner radius 

XU - upper limit on omega integration 

XL - lower limit on omega integration 

A5 - (A 1/A) **3 

A6 - 2* (C/Al) 

Y - omega 
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Parameter Study Programs 



There were three programs to complete the parameter study, and 
the program of each diffuser type is basically the same. The same 
integration scheme was used to calculate omega as the pertinent 
parameters were varied. PARSTUDY 1 was for annular diffusers, 
PARSTUDY2 for two dimensional and PARSTUDY3 for conical diffusers. 
The main device used was a set of nested do loops and the variables 
in these programs were identical to those in the previous programs 
except as follows . 

PARSTUDY1 

ALFAH - inner wall angle 

ALFAT - outer wall angle 

LDELR - dimensionless length 

TAT - tangent of outer wall angle 

TAH - tangent of inner wall angle 

DELiR - difference between inner and outer in.lpt radius 
PARSTUDY2 

LW1 - non dimensional length 
B1 - inlet dimension 
W1 - inlet dimension 

PARSTUDY3 

R - inlet radius 

LR - non dimensional length 

Program DATA was the most extensively used program in the 
and allowed the determination of all necessary numerical 
parameters for the annular diffuser data of Sovran and Klomp. 

The numerical integration scheme was used in this program to 
calculate omega, and the same symbols in the external function 
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were used as in DESIGN. The data was read from data cards and the 



parameters listed in the table of symbols were calculated. 

List of Symbols 

FAM - diffuser wall angle combination (family of Sovran and Klomp) 

AT - outer wall angle 

AH - inner wall angle 

RRAT - inlet radius ratio 

DIML - dimensionless length 

AR - area ratio 

CPI - ideal recovery coefficient 

CP - pressure recovery coefficient 

I - logical control variable 

OCF - JL c f 

RT1 - inlet outer radius 

RH1 - inlet inner radius 

DADX - area gradient 

BARL - diffuser length 

Y - omega 

RM - mean radius 

A1 - inlet area 
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DESIGN 



1 

? 
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c 



c 

•x 



4 

1 1 



6 
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EXTERNAL fct 
REAL’’- 4 L 
1 = 0 

RT] =3. 127 
P H 1 = 1 .875 
A 2=4° . f; 

AHR=1 5 . / 5 7 . 3 
A T = 1 6 . 

A TR=AT/ 57 . 3 

T A 1 P -• S I N ( ATP ) / C0S( ATR ) 

TAHR=SI M ( AMR ) / COS ( AMP ) 

C=TATR/TAHP 
n=RT 1 — R H 1 * C 
C2=D**2-A2/3.141 
B=2*P*C 
4= C *--■?- 1. 

B3=n**?-4.*A*C 2 
RH 2=1-3+83 ’• 5 ) / ( 2 . ) 

PIEL= ( RH2-RH1 ) / T A H R 
RT2='A?/' , .l4l + R H 2 * * 2 ) * * • 5 
R=DIR./CnS( ATR ) 

BIGL = ( R + DIFL )/?. 

START NUMERICAL. INTEGRATION ROUTINE 
X U = R 1 G L 
X L = 0 

A = .5< : (XU+XL) 

0=XU- XL 
r = . 484 r 8'>) *n 

V=. 0 4 G 6 3 < 1 G* ( F CT i A + C» I /\Tk » I ATR ) + I-CI l A-C , I A i k, iflHKJi 
C=. 4inui56*P 

Y=Y + . C«i 13? 4,' 3* (FCT ( A+C » T ATR , T AI1R ) +FCT ( A-C t T ATR, TAHR ) ) 

C = . 3' 668*7* D 

Y = Y + . 1 3 0 3 L 5 3 s ' : ( FCT ( A + C, TATR, T AHR ) +FC1 ( A-C , T ATR, TAHR ) ) 

C =. 1621267*D 

Y= Y + « 1861 735*{ FCTI A<C , TATR, TAHR ) +FCT ( A-C , T ATR , TAHR ) ) 
Y=p*( Y+. 1651 197* TC r( A, TATR-, TAHR ) ) 

END INTEGRATION 

W R I T L ( 6 , 3 ) Y , A T , o T 2 , R H ? , D I F L 

FORMAT ( * • , • O'-if GA= • , F * • 2 » 2 X , ‘ AL PHA= * , F5. 2 , 2X, * RT2 = ’ ,84 

I . ?,2X, • R H 2 = * , F 4 . 2 , 2 X , * D I E L = * ,F5.2) 

AT= AT+. 25 

I r { AT * GT.22 . 3 ) GU TO 4 
GO TO 2 
1=1+1 . 

WR I I E ( 6 •, 1 1 ) 

FORMAT (///, 2i -X ,' NF N AREA RATIO*, ///)■ 

I F{ I .CO. 1 . ) GO TO 6 
IF I I .E0.2. I GO TO 5 
IF (I.EQ.3) GO TO 7 
A2=5° . ( 

GO TO } 

A 2 = 6 R . 5 
GO TO 1 
STOP 
END 



FUNCTION FCTIX , TATR, TAHR ) 

RT1=3.123 A5=( XI/ ( X2-X3) )*v7. 

RHI =1 .875 A4=?.* ( { X2**. * ) + ( X3 

X l=PTl **2. F C T = A 5 * A 6 

X 2 = ( P T 1 + X * T A T 3 ) v*?, RFTU-N 

X 3 = ( P H 1 +X».« T AM? ) * * 2 , [ N ) 

//GO.SYSIN DO 



5) ) /XI 
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C CONICAL PARAMETER STUDY 

C PARSTUDY1 

// 

// 

// 



IMPLICIT REAL*4(L) 

EXTERNAL F C T 
'WRITE (6,17) 

17 FORMAT ( ' 1 • ) 

A 1 = 1 . 

R= ( 1 . / 3 . 1 A 1 ) >' # . 5 
A L P A = ]. . 

L R= 1 . 

DOl [=1,15 

DO 2 J = 1 , 1 5 

ALf AR= ALFA/57, 3 

TA\'=S IM ALF AR ) /COS (ALTAR) 

XU = L R-R 

AR=3« 141*< R + YlJ-cTAtM) **2 

CPI=l.-i./AP,**2 
X L = 0 , ( 

A = . 5* ( XU+.XL ) 

0= XU- XL 
C = . 4R4 1 BCT^n 

Y=. 0606371 9* ( FC T ( A + C, T AN ) + FCT ( A-C , TAN ) ) 

C = • 6 1 8 r 186*D 

Y = Y + . 1 3 0 3 '' 5 3 * ( FCT (A+C, TAN) + FCT ( A-C, TAN ) ) 

C=. 3R 66 867*0 

Y = V . <,<+:. 32^R .Tv ( TC I ( A + C , I AN i + FC | i A-C , l AN) ) 

C=. 162 1267*0 

Y = Y + „ 1 6 617 3 6 - ( F C T ( A+C , TAN) ♦ FCT< A-C, TAN) ) 

Y = 0*{ Y4 . 1 6 6 1 1 R 7 - T C T ( A , TAN) ) 

W RITE (6,3) Y , LR , A L F A , R 

3 FORf"\T( 1 0 y. f G A = ' , C 5.2,2X, 'L/R=' ,F5.2, 2X, • A L F A = ' ,F5.2, 

1 X, ' P=* , F 4 , ? ) 

alfa=alfa + i 
2 CONTINUE 
A L F A = 1 . 0 
L R = L 0 + 1 .0 
] CONTINUE 

S TOP 
END 



FUNCTION F C T ( X , T A N ) 

A 1= 1 . 

R= ( 1 . / 3 . 1 4 1 ) *# . 5 

A= 3. ] 41 * (R+X*T \N ) **2 

X 1= ( A 1 / A 1**3 

C=2.*3.141X(R+X*TAN> 

X 2 = C / A 1 

F C T = X 1 X 2 

RETURN 

END 

//GO. SYS IN DO * 
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c 



TWO DIMENSIONAL PARAMETER STUDY 

parstudY? 



IMPLICIT RFAL'v4(L) 

EXTERNAL FCT 
Al.FA = ( . u 
A 1 = 1 . 0 
L W 1 = ] . D 
PI=l. 

DO L M= 1 , 7 
00 2 N= ] » 7 

no 3 i=i,6 

A L F A R = A L F A / 6 7 . 1 
T AN=S I \ ( \LFAR)/COS ( ALTAR) 

W 1 = 1 . / 6 1 
L =• L W 1 * 1 

AR = Al + ?.*L*TANnU 
CP 1 = 1 . - 1 ./ARY* 2 
X U = L 
XL=C. C 

A = . 6 * ( XUf XI ) 

o=xu-xi 

C = . AB4 DV - -1*D 

Y = . i. AC 63 7 1 9* ( F f.T ( A + C , T AN f B 1 » W 3 )+FCT(A-C ,TAN,°i,Wl ) I 
C= . A 1 B r ' 1 5 6 D 

Y=Y + . 1 3030 6 ?* ( FCTI A+C , T A N » B 1 ,'.11 I+FCT ( A-C ,T AN, B1 , Wl ) ) 
C=. 31 66 857*0 

Y=Y + .t)Oi.;32^o3* ( f'CT ( A+C ,TAN, 01 , W] ) +FCT ( A-C , TAN, Bl , Wl )l 
C = . 1.62’ 267*0 

Y= Y + . ] Sol 7 3 6' ( «CT ( A+C , TAN, Bl , W] ) +FCT ( A-C , T AN , B 1 , W 1 ) ) 

Y ( y . i h 6 i l y 7* f-C I ( A , i AN ,ui, W i ) ) 

WRIT! { 6 ,4) Y, LW 1 , Wl , PI , ALT A 

6 FORMAT ( * 0‘1 EG A = • , F6.2, ?X, • L/W l=',P5.2,2X,*Wl=*,F5.2,2X 

1 , 'Bl=« , 1*4. 2 ,2X , ' ALP A = « , FA. 1 ) 

B 1 = A 1 +• I . 

3 CONTINUE 

B 1 = 1 . 

ALFA=AI FA *2 .6 
2 CONTINUE 

B 1 = 1 . 

ALT A= 1 „'■) 

L Wl = l.Wl *2. 

I CONTINUE 

STOP 
END 



FUNCTION FCTIX , TAN,B1,K1 ) 
A 1 = 1 . 

A = A 1 + ? . * X y T A N Y B 1 

X ] = ( A l / A ) y * 3 

C = 2 . Y w 1 + 2 . Y B 1 + 6 . * X y T A N 

X 2 = C / A 1 

FCT=X1 YX2 

R n TURN 

END 

//GO. SYS IN 00 * 



74 



c 



ANNULAR PARAMETER STUDY 
PARS TIJD Y3 



C 



IMPLICIT R E A L * 4 ( .L ) 

EXTERNAL f-CT 
WRITE ( 6 , 2*> ) 

2 5 FORMAT! *1') 

A 1 = 2 . i . 

ALFAH = <'.D 

ALF AT= 1 . 1/57.3 

L I 7£LR= 1 .0 

R I = 2 . 0 

Of] 1 M=l,7 

DO 2 J = 1 ,21 ,5 

N = J + 1 

DO 3 K = N , 2 2 , 5 
DO 4 I=],4 

R0= ( A 1/ 3. 14 1+P I **2 ) **. 5 

DhLR=PO-RI 

A I E A7R = ALE A T*5 7 . 3 

At FAUR= ALE AH 45 7 . 7 

TAOS IN! ALE AT ) /CDS ( ALE AT ) 

1 AH=S Ifl ( A l F AM) /CDS < ALF AM) 

XlJ = l DL L R *DE L P 

A P = ( 3.141*! ( RP+XU*T AT ) 4*2- ( R I +XU*TAH ) **2 ) ) /A I 

C P l = l . - ( 1 . / A R ) * * 2 

PA r >X = 2. 4RU+TAT -?.*Rl*TAH 

D 1 F F " ( ( TAT ) **2 - ( TAM )**2 ) *2. 

XL =0.0 

A- . 5'" i XU-t-XL ) 

OXU-Xl 

C= . 4 0 4 ' 50] *D 

Y = . 0 4 ’• 6 3 7 1 0 * ( F C T ( A + C , R I , RO , T A T , T AH ) + FC 1 ( A-C , P I , RO , T A T , 
1 TAM ) ) 

0.4] 8" 1 56*0 

Y = Y+ .''0324' P* ( FCT ( A + C , R I , RO, T AT , TAM) + FCT! A-C ,R I , PO, TA 
IT, TAM)) 

C = . ?■ 64 857*0 

Y = Y * . 1 3 M 3 f ■ 5 3 * ( t-'CT ( A+C ,RI , Ru , T AT , T AH ) +F C T ( A-C , R I ,R0, TAT 
1 , TAM) ) 

C = . ) 6? 1 2 67* n 

Y -Y + . 1541 7 3 5 * ( F CT ( A +C , R I . RU t T AT , T AH ) + FCT ( A-C * R I , R 0, T A T 
1 , TAM) ) 

Y=U* ( V + . ] Of. 1 1 9 7*FC J ( A , R I , RO , T AT , T AH ) ) 

KRITF (6 , 8) Y ,LDf LR, DADX 

5 FUR MAT (• « , *U’! F: C, A = ' , F5 . 2 , 3 X , • L D F LR= • , F 6 . 2 t 3 X , *DADX= ' , 

1 F 5 . 2 ) 
k I = R I + ? . 0 
4 CONTINUE 

R 1 - ? • 0 

ALFAT=ALFAT+5./57. 3 

3 CONTINUE 



2 



c* 0 2.0 

ALFAH = ALOH+5. / 57. 3 
At FAT =ALEAH + 1 . / 57. 3 
CONTINUE 
R I = 2 . 0 



FUNCTION f : C T ( X ,RI ,FO,TAT , T A hi ) 

A 1 = 2. . 

X2 ~ ( RO+X*T AT )*■ ? 

X 1 = { P I -# X*T\F)* 1 2 
A 7 = 3 . 1MMXI-X] ) 

X 3 = ( M **2)/ (A3** 7 ) 

X4=2.*3. 141*! ( XI*-*. 5) + ( X2**. 5) ) 
F C T = >' 3 * X 4 

RETURN 

END 



1 



Al. F AH = C .0 

A L F A T = 1 .0/57.3 

LDELR=LDtl R +2. 

CONTINUE 

SI OP 

END 



.SYS IN o-i 
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DATA P. EDUCT 1 0 s ! FRO M SOVRAN ANO KLOMP 



EXTERNAL FCT 
PFM*4 l 
RT1 = 7. 55 
WRITE (6,45 ) 

45 FORMAT ( ' 1 ' ) 

1 R F A D ( 5 , 10) FAM , AT, AH, RRAT, Dl Ml. , AR,CPI ,CP, I 

1C FORMAT (3 F5. 2, F 3. 2, F 5. 2 , F 5 . 3 , 2F 4 . 3 , I 1 ) 

OC p = C.PJ-CP 
PHl = RTl*OR. AT 
DELR=RT l-RHl 
A TR = A T / 5 7* 3 
AHR=AH/ 57. 3 

T AT F ~ S 1 N ( \TR I / COS ( A TR ) 

TAliR-S I r i { A H R ) / COS ( A HP ) 

UAl)X=?. * (R r 1 ,; = T A T R- P H 1 * T At I R ) 

BARL=DIML* { RT l-RHl ) 

XU=RARl 

XL=9«C- 

A = . 5 * ( x U * X L ) 

D=XU-XL 
C = . 4 8 4 ‘ 3 c 1 * n 

Y= . f'4'.i 427 1 °- ( F CT { A+C, ATR, AHR , RH1 ) f FCT ( A-C , ATR, AHR , RH] 
1 ) ) 

C=.41 3 015 6 * D 

Y = Y +- . i 9 0 3 2 40 8 * ( F C T ( A+ C , AT P , A HR , K H 1 ) + F C T ( A- C , A 1 R , AH R , R 
1H1N 

C-. 3 u 66857* n 

Y=Y+. 1 3F3-' 53* ( FCT < A+C . ATR . AHR. RHD + Ff T ( A-C . ATR . AHR . PH 
ID) 

C = . 5. 62 l 26 7*n 

Y = Y+. 1 5M 7 3 5*( FCT ( A+C , ATP , AHR, PHD ) 

Y = D - ' ( Y+ . 1 651 1 0 7vfCT( A , ATR , A HR, RH1 ) ) 

PARI- { 0 A n X - l . 5 ) / DEL R 

PARI 1 = .,CF/ PARI 
PAR2=( DADX/DtL R ) *'2. 

PAR2DHCF/PAR2 

P A R 4 A = 2 4 ( T ATO-RK AT * TAHR ) / ( 1 . +RR AT ) 

P AR 5“ ( C R/C D I ) yp AC 4A 
CF = r>AR4 A* ( D-CP/CP D 
WRITE { 6 , AO) r A ‘ i , P A R 4 A , C F , P A R 5 , Y 
90 FORM A T ( * « , TA';= ' , F0.2, 3X, »PAP4A=' ,L r l'.3,3X, 'CF=* ,510 

l . 3 , 3 X , ' C P / ' ) v F G A = ' , t 1 0 .3,3 X , 'UMEGA=» ,F5.2) 

V. R I IT ( 4,77) 

77 FORMAT (/) 

I F ( I .ED. 1) GO TO 4 
GU TO 6 

4 WRITE (6, 7) 

7 FORMAT (///, ?f. X ,' MEW FAMILY',/-//) 

6 I F( I. FG.2) GO TO 8 

GO TO 1 
3 STOP 

END 



FUNCTION F C T { X , ATR , A FR , RH 1 ) 
R T 1 = 7 . 5 5 

X 1 = P T 1 *=:2-PHl**2. 
f A T P = S 1 f ' ( ATT 1 /C 0$ ( ATR ) 

T AHR= S I M ( AH P I /COS ( AHR ) 

X? = { FT OKO ATR ) •' - ? . 

X 3= ( P 0 1 +• X ' c T \M0 ) **? . 

A6=( Xi/ ( X2.-X3) ) 3. 
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A6 = ?« * ( ( X2#*. 5 ) + ( X3**. 5 ) ) /XI 

FCT=A5*A6 

UPTURN 

F:ND 

//GO. SYS IN DD * 
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